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Spin and Valley Splittings in Multilayered Massless Dirac Fermion System 
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The inter-layer magneto-resistance in a multilayered massless Dirac fermion system, a-(BEDT- 
TTF)2l3, under hydrostatic pressure was investigated. We succeeded in detecting the zero-mode 
(n=0) Landau level and its spin splitting in the magnetic field normal to the 2D plane. We demon- 
strated that the effective Coulomb interaction in the magnetic field intensifies the spin splitting of 
zero-mode Landau carriers. At temperatures below 2K, magnetic fields above several Tesla break 
the twofold valley degeneracy. 
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Since Novoselov et al. [1] and Zhang et al. Q 
experimentally demonstrated that graphene is a zero- 
gap system with massless Dirac particles, such systems 
have fascinated physicists as a source of exotic systems 
and/or new physics. At the same time, the quasi- 
two-dimensional (Q2D) organic conductor ev-(BEDT- 
TTF) 2 l3_ (BEDT-TTF=bis(ethylenedithio) tetrathiaful- 
valene) Q was found to be a new type of massless Dirac 
fermion (MDF) system under high pressure [J-Q- In con- 
trast to graphene, this is the first bulk (multilayered) 
zero-gap material with Dirac cone type energy disper- 
sion. Another important difference is the fact that ac- 
cording to the band calculation [4, 5] or the first princi- 
ples band calculation [6j, the present system has tilted 
Dirac cones. Hence, the Fermi velocity, i>f, is highly 
anisotropic. Thus, a-(BEDT-TTF)2l3 provides us a new 
type of MDF system with a layered structure and an 
anisotropic vp. 

One of the characteristic features of transport in the 
multilayered MDF system is clearly observed in the inter- 
layer magnetoresistance shown in Fig. 1(a) j9|. The re- 
sults are interpreted as follows. 

In the magnetic field, the energy of Landau lev- 
els (LLs) in zero-gap systems is expressed as -E n LL = 
±y / 2e?it;p|n||i?|, where n is the Landau index and B is 
the magnetic field strength [10]. One important differ- 
ence between zero-gap conductors and conventional con- 
ductors is the appearance of a (n=0) LL at zero energy 
when magnetic fields are applied normal to the 2D plane. 
This special LL is called the zero-mode. Since the energy 
of this level is Ep irrespective of the field strength, the 
Fermi distribution function is always 1/2. It means that 
half of the Landau states in the zero mode are occupied. 
Note that in each LL, there are states whose density is 
proportional to B. The magnetic field, thus, creates mo- 
bile carriers. 

For k B T < -Eill, most of the mobile carriers are in the 
zero-mode. Such a situation is called the quantum limit. 
The carrier density per valley and per spin direction in 
the quantum limit is given by D{B) = B/2(f>o, where 
<f>Q = h/e is the quantum flux. The factor 1/2 is the Fermi 



distribution function at Ep . This effect is detected in the 
inter- layer resistance, R zz , in the longitudinal magnetic 
field. In this field configuration, the interaction between 
the electrical current and the magnetic field is weak be- 
cause they are parallel to each other. Hence, the effect 
of the magnetic field appears only through the change in 
the carrier density. In this regard, the large change in 
the density of zero-mode carriers gives rise to remarkable 
negative magnetoresistance in the low magnetic field re- 
gion, as shown in Fig. 1(a). It is written as R zz oc 1/B. 
This behavior agrees well with the realistic theory of Os- 
ada [Il|. 

At a high magnetic field, on the other hand, the Zee- 
man effect plays an important role in zero-mode carriers. 
In a magnetic field, each LL is split into two levels with 
energies -E n LL ± A_B s /2, where AE s =g* /i^B is the spin- 
splitting and g* is the effective g-factor. This change in 
the energy structure gives rise to a change in the car- 
rier density in LLs. At a low temperature and/or a high 
magnetic field, where k B T < AE S , this effect becomes 
important because the energy level is shifted from the 
position of Ep. It works to reduce the density. The mag- 
netoresistance changes from negative to positive, and is 
dependent on the magnetic field and the temperature as 



R zz ex 1/D(B) ■ exp( g y B B/2k B T). 



(1) 



Hence, R zz ■ B obeys the exponential law as shown in Fig. 
1(b) at a high magnetic field. 

The detection of the zero-mode and its spin-split level 
at a low field strongly indicates the high purity of our 
system and leads our investigation to a new stage. Each 
zero-mode in this system has twofold degeneracy (valley 
degeneracy) originating from two Dirac cones in the first 
Brillouin zone [4-6]. According to the simple Storner- 
like theory of quantum Hall ferromagnetism [12j . the 
Coulomb interaction plays an important role in spin and 
valley symmetries in the magnetic field. In this paper, we 
experimentally demonstrate that the valley degeneracy 
in this system is broken at magnetic fields above several 
Tesla at temperatures below 2 K. In graphene, on the 
other hand, this state was realized in the magnetic field 



above 20 T [13j. This critical field strength depends on 
the scattering broadening energy of the zero-mode. The 
broadening energy of the LLs in our system is much lower 
than that in graphene. Another important fact is that 
according to the theory of Kobayashi et a/., when the 
Dirac cone is tilted, the long-range Coulomb interaction 
has finite matrix elements between valleys, and it breaks 
the valley degeneracy in the magnetic field (Fig. 2) jl4| . 

We succeeded in obtaining evidence of the valley-split 
levels in our material, as follows. 

In the magnetic field, the effective Coulomb interaction 
enhances the spin splitting as, 



B(T) 



AE S = gf j, B B 



(2) 



where g ~ 2 is the g-factor, / ~ e 2 /eln ~ 50vB/e meV 
is the e ffectiv e Coulomb interaction in the magnetic field, 
In = y/h/eB is the magnetic length, and e is the effective 
polarizability. Then, we have the effective g-factor, 



9 



AE s /fi B B = g 



(3) 



It should exceed 2 and depends on the magnetic field 
because I ex y/~B. 

This effect can be seen in R zz under a longitudinal 
magnetic field (Fig. 1(a)). At the minimum, we have 
g* \i-qB ~ r, where T is the scattering broadening en- 
ergy or the thermal broadening energy of each LL [111 ]. 
At a high magnetic field, on the other hand, R zz obeys 
an exponential law as Eq. ([TJ. Note that R zz can de- 
tect only the charge gap, AE, between levels that are 
nearest to Ep, assuming that E-p always locates at the 
contact point. In the valley degenerate state, hence, 
AE = AE S = <?*/ib-B. The magnetoresistance in the 
valley splitting state, however, is characterized by the 
charge gap between (1, — ) and (f, +), AE = AE S - A v , 
where A v is the valley splitting, and we describe the level 
structure of the zero mode using spin index (f, 1) and val- 
ley index (+,—), as shown in Fig. 2. 

In our previous work, we published data for R zz at 
magnetic fields of up to 7 T and succeeded in detecting 
the zero- mode and its spin-split level |9| . To demonstrate 
that the magnetic fields break the valley degeneracy in 
this system at low temperatures, however, a higher mag- 
netic field becomes important. In this work, we examined 
the magnetoresistance at magnetic fields of up to 14 T in 
the temperature region below 4.1 K. 

A sample to which four electrical leads were attached 
was encased in a Teflon capsule filled with pressure 
medium (Idemitsu DN-oil 7373). The capsule was set 
in a clamp-type pressure cell made of MP35N hard alloy 
and hydrostatic pressure of up to 1.7 GPa was applied. 
Resistance measurements were carried out using the con- 
ventional DC method with the electrical current applied 
along the c-crystal axis, which is normal to the 2D plane. 

Figure 1(a) shows the magnetic field dependence of 
R zz below 4.1 K. It reproduces well our previous data 
published in Ref. Q. We can understand this magne- 
toresistance in terms of zero-mode carriers, including the 
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FIG. 1: (a) Magnetic field dependence of R zz below 4.1 K 
under pressure of approximately 1.7 GPa. (b) R zz in the 
low field region below 0.4 T. (c) Magnetic field dependence of 
R zz ■ B. (d) Temperature dependence of B p (solid triangles) 
and Bmin (solid circles). We show B p and Bmin estimated 
from previous data [3] by open triangles and open circles. 
The solid line is the curve of Sill with df~4x 10 4 m/s. 
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FIG. 2: (color online). Zero- mode and its spin- and valley- 
split levels for the cases AE S > A v (a) and AE S < A v (b). 
The level structure is described using spin index (f, J.) and 
valley index (+, — ) 



spin splitting. In the interpretation of R zz in Ref. [9|], 
however, we did not take into consideration the Coulomb 
interaction. The Coulomb interaction should be taken 
into consideration because it has a significant influence 
on the transport phenomena in the magnetic field. 

Thus, the first step is to verify the effective Coulomb 
interaction in this system from the effective g-factor g^- ln 
estimated from the relationship 3* nin /iB-B m i n ~ T at the 
magnetoresistance minimum. Here, broadening energy T 
should be proportional to k B T because -B m i n decreases 
linearly with decreasing temperature, as shown in Fig. 
1(d). It can be written approximately as T = To + 2k B T 
from the investigation of the relationship between T and 
the temperature from the simulation calculation of R zz 
along Ref. [ll|. The scattering broadening energy, r , is 
roughly estimated to be about 3 K as follows. 

Since each LL is broadened by the scattering of carriers 
and/or thermal energy, the zero-mode is sure to overlap 
with the other LLs at a low magnetic field. In such a 
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FIG. 3: (color online), (a) Temperature dependence of g^ in 
(solid circles) and g* xp (squares), (b) Magnetic field depen- 
dence of AE B and g^ n . Each data can be reproduced well 
with Eq. © (blue line) or © (red line), with e ~ 190. 



FIG. 4: (color online), (a) Temperature dependence of AE 
(squares) and A v (triangles) for B= 7 T. Solid line is AE S ~ 
18 K prospected at 7 T. (b) AE for B=7T below 1 K. 



region, the relationship of R zz ex 1/ B loses its validity. 
We can recognize this region in Fig. 1(b) where a pos- 
itive magnetoresistance is observed. This critical field, 
B p , shifts to a lower field with decreasing temperature 
down to about 2 K, where it almost saturates at about 
0.04 T, as shown in Figs. 1(b) and (d). The overlap be- 
tween the zero-mode and other LLs, primarily the n=l 
LL will be sufficiently small above B p and as a result, 
the negative magnetoresistance is observed there. Then, 
we have a tentative relationship: -Eill ~ ^k^T v at B p . 
In fact, Sill with np' v 4x 10 4 cm/s is reproduced well 
except in the temperature region below 2 K. This Fermi 
velocity vf corresponds to that estimated from the tem- 
perature dependence of the carrier density, N, written as 
N ex (T/v F ) 2 , with v F ~ 10 5 m/s within a factor of 3 0. 
The discrepancy of the data from the curve of Sill below 

2 K, on the other hand, suggests that thermal energy is 
sufficiently lower than the scattering broadening energy 
To. Thus, To is roughly estimated to be app roximately 

3 K from the constant value of B p as Tq = y/2ehvp\B p \ 
[15j . Note that in the estimation of Tq, we did not con- 
sider the spin splitting because r has little impact. This 
scattering broadening energy is much lower than that of 
graphene. In graphene, Tq was estimated to be about 30 
K 0. 

Using the relationship g^ in {i B B m i B ~ T, g^ in at each 
temperature is estimated as shown in Figs. 3(a) and (b). 
It gives us definite evidence that the effective Coulomb 
interaction intensifies the spin splitting as follows. 

First, at each temperature or magnetic field, g^ in ex- 
ceeds 2. It increases by about three times from about 
4.5 at 4.1 K to nearly 12 at 0.1 K, as shown in Fig. 
3(a). Note that this strong temperature dependence of 
9min ma y arise from the effect of its strong magnetic field 
dependence because B m - m depends strongly on temper- 
ature, as shown in Fig. 1(d). In fact, we can find this 
effect in the magnetic field dependence of g^- ln shown in 
Fig. 3(b). In the case of a system without the effective 
Coulomb interaction (1=0), g* should always be 2 in the 
magnetic field. However, g^in depends strongly on the 
magnetic field, as shown in Fig. 3(b). It obeys 1/y/B. 
Assuming that e is independent of temperature and using 
Eq. ©, with e ~ 190, we try to fit the curve in Fig. 3(b). 



This simple formula reproduces the data well and shows 
evidence that the effective Coulomb interaction plays an 
important role in the spin splitting of this system. In 
particular, at a low magnetic field, this effect is strength- 
ened. Here, we note that e ~ 190 is the effective polar- 
izability when the Coulomb interaction engages directly 
in the spin splitting. Thus, this is the upper limit of this 
system. Recently, Morinari and Tohyama performed a 
simulation calculation for in-plane magnetoresistance us- 
ing this effective polarizability [17|, and our anomalous 
data [7[ are quantitatively reproduced. 

To conclude, we can detect the effect of Coulomb inter- 
action probed by inter-layer magnetoresistance measure- 
ments in this system. According to the simple Storner- 
like theory of quantum Hall ferromagnetism [12j . the 
strong Coulomb interaction leads us to summarize that 
the twofold valley degeneracy may be broken at a high 
magnetic field that satisfies y/Tr/21 > T. 

The last step is to determine the effect of valley split- 
ting in the magnetic field. At 7 T, for example, the ef- 
fective Coulomb energy exceeds T below 2 K. Thus, we 
expect that this effect appears in the effective g-factor, 
g*, estimated from the resistance curve obeying Eq. ([T]) 
as shown in Fig. 1(b), at temperatures below 2 K and 7 
T. 

At temperatures above 2 K, g* ~ g^ini as shown in 
Fig. 3(a). The slight difference arise from the effective 
Coulomb interaction obeying Eq. ([3]). At 7 T, for exam- 
pie, AE = g* xp ^ B B ~ AE S ~ 18 K (Fig. 4). 

The behavior of g* and AE below 2 K is much more 
impressive. Both abruptly drop by approximately 1/50 
times from g* xp ~ 3.8 and AE ~ 18 K at 2 K to g* xp ~ 
0.08 and AE - 0.4 K at 0.07 K. Judging from the detec- 
tion of a strong effective Coulomb interaction, we ascribe 
the drop of g* and AE in Figs. 3(a) and 4 to the break- 
ing valley degeneracy in the magnetic field above several 
Tesla, as shown in Fig. 2(a). In graphene, it was real- 
ized at magnetic fields above 20 T [12| and the realistic 
theory established that the strong Coulomb interaction 
broke the valley symmetry at a high magnetic field [13| . 
In our system, on the other hand, T is much lower than 
that of graphene and therefore, we could detect the effect 
of valley splitting at a low magnetic field. At 2 K, for ex- 



a mple , the critical magnetic field strength that satisfies 
y/t/2I > T is about 7 T. At f K, it is about 4 T. 

Based on this speculation, valley splitting in our sys- 
tem can be estimated tentatively as A v = AE S — AE. 
At 7 T, for example, valley splitting abruptly emerges at 
about 2 K and increases up to approximately 17 K at 
0.07 K, as shown in Fig. 4(a). 

Recently, Kobayashi et al. predicted the pseudo-spin 
polarized ferromagnetic state with an easy plane where 
A v > AE S as shown in Fig. 2(b) 14]. In this state, 
excited vortices and anti-vortices can move and the sys- 
tem may undergo the Kosterlitz-Thouless (KT) transi- 
tion to form pairs of vortices and anti-vortices at the 
lowest temperature. Our data, however, suggest the sit- 
uation A v < AE S , as shown in Fig. 2(a). 

Examining AE data below 1 K in more detail (Fig. 
4(b)), we find that it seems to go across zero and then 
the pseudo-spin polarized ferromagnetic state in which 
the levels between (i, — ) and (f, +) cross may be realized 
below 0.07 K (Fig. 2(b)). We also expect that this state 
would be realized at a high magnetic field. The data be- 
low 0.8 K in Figs. 1(a) and (b) are apparently discrepant 
from Eq. ([I]) at a high magnetic field. The saturation 
of R zz may indicate a symptom of the crossover to the 
pseudo-spin polarized ferromagnetic state at a high mag- 
netic field, as shown in Fig. 2(b). 

Lastly, we briefly mention the mobility edge in this sys- 
tem at a low temperature and a high magnetic field. In 
Fig. 4(b), we perceive the fact that AE below 0.1 K is 
lower than To ^ 3 K and yet, R zz at 7 T is much higher 
than that at the minimum. At 0.07 K, for example, AE 
is approximately 0.4 K. It arises from the localization of 
electrons on the tails of the broadened zero-mode. When 
LL is broadened by scattering, only the electrons at the 
vicinity of center i? n LL are mobile. This critical energy 
_E n LL ± E c is called mobility edge. The smaller the Lan- 
dau radius is, the stronger the localization is. The effect 
of localization also depends on the temperature. Assum- 
ing that the mobility edge in this system locates at the 



vicinity of the center of each Landau level, E c should be 
less than AE/2 ~ 0.2 K at 0.07 K. 

In conclusion, we succeeded in detecting the zero- 
mode and its spin- and valley-split levels in the multilay- 
ered MDF system a-(BEDT-TTF) 2 I 3 probed by inter- 
layer magnetoresistance measurements. The effective 
Coulomb interaction plays an important role in inten- 
sifying the spin splitting of zero-mode carriers. Hence, 
the effective g-factor exceeds 2 and depends strongly on 
the magnetic field as g* ex \j\[B. The twofold valley 
degeneracy, on the other hand, is broken in the magnetic 
field that satisfies yJir/21 > T. At 7 T, valley splitting 
abruptly emerges at about 2 K and increases up to ap- 
proximately 17 K at 0.07 K. This system is very pure and 
therefore, we could detect the effects of Coulomb inter- 
action and its associated valley splitting in the magnetic 
field, which is much lower than that of graphene. 

Recently, the multilayered MDF system was also re- 
alized in one of the parent compounds of the oxyp- 
nictide superconductors [18|]. This MDF system, how- 
ever, is not simple because this compound has multi- 
band structure and therefore some Fermi surfaces. In 
this sense, a-(BEDT-TTF)2l3 provides a suitable test- 
ing ground for the transport of the multilayered MDF 
system. According to the theory by Shon and Ando, 
the mono-layered MDF system show a universal conduc- 
tance, a xx ~ 2e 2 /ir 2 h, which is independent of the mag- 
netic field, the temperature and the scattering at Dirac 
point [lj|. In the multilayered MDF system, however, the 
tunneling between layers reflects the broadening of zero- 
mode and it loses the validity of universal conductance in 
the magnetic field. Realistic theory reproduced well our 
anomalous in-plane magnetoresistance [3, H7j| . Further 
investigation will lead us to new phenomena character- 
ized by the multilayered MDF system. 
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